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fvq I Abstract The electron density n{r,t), which is the central tool of time-dependent density 

functional theory, is presently considered to be derivable from a one-body time-dependent 
potential V(r,r), via one-electron wave functions satisfying a time- dependent Schrodinger 
equation. This is here related via a generalized equation of motion to a Dirac density matrix 
now involving t. Linear response theory is then surveyed, with a special emphasis on the 

^ , question of causality with respect to the density dependence of the potential. Extraction of 

j^ ■ ^(fiO for solvable models is also proposed. 
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^1 ■ 1 Background 

Early work on time-dependent density functional theory can be traced back at least to 1972 
ni l2][3t followed by various studies I4ll5] |6ll7ll in the 1970s / early 1980s, and culminating in 
the important proposal of Runge and Gross (RG) fSl. In essence, the RG argument general- 
^Q ' izes the Hohenberg-Kohn theorem |9| to time-dependent external potentials. Though parts 

f~^ , of this important study were questioned in 1 101, the RG work is widely accepted as the basis 

for the assertion that, for a specified initial state, there is a unique correspondence between 
the time dependent density n(r,f) and the external potential Vext(r,^). This points the way to 
construct a time-dependent extension V(r,f) of a static Slater- Kohn-Sham (SKS) like poten- 
5h ' tial y(r) 111 . 12 1 which will then generate a Slater determinant of non-interacting electron 

wave functions, (^/(r,f) say, satisfying the time-dependent Schrodinger equation 

~V^ + V{r,t)\Ur,t) = ihj^^i{r,t), (1) 
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with a specified determinant at time t = 0. Then the electron density n(r, t) introduced above 
is constructed, formally exactly for N electrons, as 

N 

n{r,t)=Y.^;{r,t)^i{r,t). (2) 

Of course, formal exactitude requires precise knowledge of the one-body potential V{r,t) 
in eqnITI. At the time of writing, such knowledge is limited for the key exchange and corre- 
lation contributions entering V (see also eqn([9]l below). These can be formally constructed 
[13"14','151, and perturbative approaches that converge on their exact form are also known 
06, 17, 18|. In practical terms, however, accurate exchange-correlation functionals beyond 
the local density approximation in space and time are still elusive and remain a topic of 
general interest. 



2 The challenge of Schirmer and Dreuw 1 10 1 to the RG arguments, and some 
responses 

This is the point to return to the work of IIIOI . This study contains serious criticism levelled 
against the very foundations of TDDFT (articles I19II20||2T|| are also concerned with the 
basis of the theory). In fTO|, the variational basis of TDDFT proposed by Runge and Gross 
(§1 was not only challenged but seemingly refuted. To be more specific, Schirmer and Dreuw 
claimed that the variational derivation of the time-dependent SKS equations in (H is not 
valid due to an ill-defined action functional proposed there. A non-variational treatment 
would also encounter difficulties, since in this case the SKS system would permit one to 
reproduce, but not predict, the exact electron density. 

Two contributions involving the present authors I23II24I have been motivated by the 
criticism in 1101 of the RG work. Both of the contributions accept the challenges of the 
RG proof, but do not require one to abandon the RG conclusion nonetheless. Let us start by 
summarizing the content of 1 23 ], because this is very specifically focused on time-dependent 
theory, whereas |24|, though also motivated by the challenges in 1101 , is basically dealing 
with time-independent DFT. 



2.1 Solvable example of a family of two-electron model atoms with general inter-fermion 
interaction: dynamical generalization 

As brief background to the above example, Holas, Howard and March (HHM) [25}, ob- 
tained analytical solutions for ground-state properties of a whole family of two-electron 
spin-compensated harmonically confined members characterized by a given interfermionic 
potential energy u{ri2) (See also I26II27II ). In L23J a start is made on the dynamic generaliza- 
tion of the harmonic external potential. In the above context, a simplified expression is ob- 
tained for the time-dependent electron density for arbitrary inter-particle correlation, which 
is completely determined by a one-dimensional non-interacting Hamiltonian. That such a 
construction is generally possible has been shown by Qian and Sahni I13II14I . but follows 
also from the harmonic potential theorem f28"291 for this specific example. Furthermore for 
the simplest case: the Moshinsky atom |30], where the interaction u{ri2) is also harmonic, a 
closed solution for the Fourier transform of the density, namely the time-dependent atomic 
scattering factor, is found. 
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To summarize the essence of the time-dependent density n(r, t) calculation in I23II . from 
the above model, we take the special but nevertheless important case of a system which is in 
its ground-state at f = 0. After generalizing the static separation of center of mass (CM) and 
relative motion (RM) to the dynamic example under consideration, the above assumption at 
t = leads to the square of the CM wave function as the simple Gaussian form 

where the time dependence is determined by the length scale acM(0 of the oscillator. 

For the Moshinsky example I30II the time dependent atomic scattering factor f{k,t), 
defined by 

/(k,f) = |«(r,Oe*'rfr, (4) 

is the convenient tool. The total scattering factor turns out, for u{ri2) = —^J(fr\2, to have 
the form 

fiotik.t) = 2/cM(A:,mcm<?(0)/cM(A:/2,wJcn,|(f))- (5) 

where acu_{t) entering equilS} is related to m^j^^it) in eqn© by 

acM(0 = ITT- (6) 



Or more generally, the dynamic generalization of the static HHM density is obtained in [23 
in terms of the time-dependent relative motion wave function as 

"("■'0 = ^^exp( 2 



«Cm(0 

:j^"rfy/exp(-^)|v/or^°(flCM(Oy,of 
sinh(ry/acM(0) 



(7) 



('7/«cm(0) 

Though this is an admittedly simplistic two-electron correlated time-dependent prob- 
lem, the time-dependent density n{r,t) can be got via a one-body time-dependent potential 
y(r,f), thereby supporting the original RG assertion. There is no conflict either, we hasten 
to add, with the Schirmer-Dreuw conclusions. These authors, in spite of questioning the RG 
derivation fSl of the SKS equations, nowhere claim to have disproved this important asser- 
tion! On the contrary, an alternative proof is provided in 131] that the mapping is indeed 
valid. 

To conclude this sub-section, we stress that a correlated two-electron example proposed 
in the static limit in 1251 has been solved exactly in the dynamic generalization in which the 
system is in its ground state at time f = 0. In particular, eqn(|7]i allows the time-dependent 
density «(r,f) of the correlated dynamical problem to be reduced to the single-particle 
problem of calculating, probably numerically, from a one-body time-dependent Schrodinger 
equation, the relative motion wave function V'ooo ^ ■ "^^^ wave function, though calculated 
from a one-body equation, involves the sum of the harmonic confinement potential and the 
interparticle potential energy u. Having obtained the exact density, also the generating SKS 
potential may be obtained by a straightforward, but numerically non-trivial, inversion of the 
SKS equations in this two-electron case |32 33]. The one-to-one correspondence of densi- 
ties and potentials in the non-interacting case can hence be verified quite explicitly for this 
specific example. 



2.2 Linear response theory and its inversion 

Since the comments |24| are basically time-independent, we shall defer these to follow a 
brief discussion of the important linear response function in time-dependent theory, already 
to the forefront in the discussion in QJ and ^. For the non-interacting SKS system gener- 
ated by potential V{r,t) with first order self-consistent change AV{r,t), with corresponding 
response function X!(r,f,r',r')' we have for the density change An{r,t) the formal result 

An(T,t) = fdt' [ dr'xs{r,t,r',t')AV{r',t'). (8) 



Here the first order change AV in the (now time-dependent) SKS potential is given by 

AV{r,t) = V,,,{r,t)+ J dr'^0^+ J dt' J dr'f,,{r,ty,t')An{r',t'), (9) 

where /^c is the as yet unknown, exchange-correlation kernel I34II35I . 

In this last cited references by Gross and coworkers, the issue of causality was first 
raised, which was also taken up later by Amusia and Shaginyan I36II37 1 and Harbola 1381 
[39', 40 1 . The question posed is quite general: namely whether the potential depends on the 
density in a causal manner. This, we believe, is another question important to the foundations 
of TDDFT under discussion here. 

The causality issue is important in its own right, but has also implications for the varia- 
tional formulation of TDDFT: The action principle proposed by Runge and Gross [8] leads 
to a symmetric and hence unexpected non-causal form of the inverse of the response func- 
tion X, giving rise to what has been termed the symmetry-causality paradox (see also [41]). 
Different reformulations of the action principle have appeared in the literature, e.g., |42l!43l. 
up to a recent contribution by Vignale in which the principle of least action in its conven- 
tional form is abandoned |^4l. 

Returning to the question of causality, we follow the treatment of Amusia and Shaginyan 
1371 who write the external potential in terms of the density, by invoking the many-body 
linear response function x (in contrast to the one-body Xs used in eqn([8]l above). With the 
external potential Vext(r,f) one can write 

An(r,t) = j dt' j dr'x{r,t,r',t')V,n{r',t'). (10) 

In I37II . it is assumed that eqn jlOl l can be changed to the Volterra integral equation (see also 

my- 



1 d^An{t) _^^ ,^, , f^, 1 d^x{t,n, 



'-^^^\h'm-^'-^''^^ 



(11) 



K{t) dt^ exiv y . y ^^^^ ^^j 

where the spatial variables are omitted for clarity of notation. The function K{t) entering 
eqnilllb is defined by 1371 

dx{t/) 



m 



dt' 

The solution of eqndllt has the form I37II40II 



(12) 



„ , , 1 d^AnU) f ,., „ 1 d^A,i(t') 

where R(t,t') vanishes for t < t'. Harbola |40| argues from the above treatment that while 
the density depends on the potential in a causal manner, the reverse is not true, leaving 



no room for a symmetry-causality paradox. On the other side, Amusia and Shaginyan 1361 
[37l as well as Cohen and Wasserman |41 ] provide arguments that causal inverse response 
functions indeed exist. Our present contention is that all-important for a decisive answer to 
the causality question resides in the mathematical properties of eqnill3l>. This is an important 
future area therefore for rigorous mathematical physics. 



2.3 Two further comments pertaining to the Schirmer-Dreuw study 

At this point, we return to the comments of Holas et al. |24| on the Schirmer-Dreuw study 
1101 . As already mentioned above, in [24] the two particular subjects tackled are concerned 
with the original, time-independent DFT but have implications also for the time-dependent 
generalization. In the present section, we merely summarize the relevant points in Ii24] for 
the present context. 

The first focus of |24| was to answer in a positive fashion a question posed in flOI . 
Whether a local operator can be reconstructed from knowledge of its particle-hole (p-h) 
matrix elements when the number of particle states exceeds one is the essence of the ques- 
tion. It arises in the context of a linear response treatment within TDDFT, in which only 
p-h matrix elements of the perturbing operator appear in the relevant equations, although an 
exact many-body approach requires also the knowledge of p-p and h-h elements. It turns out 
that there is no conflict with the exactness of TDDFT due to this apparent loss of informa- 
tion 1 10|. Schirmer and Dreuw formulate and prove the theorem that a local (multiplicative) 
operator u = u{r) is uniquely determined to within a constant by its p-h and h-p matrix el- 
ements with respect to a complete one-particle basis and any partitioning of that basis into 
occupied and unoccupied one-particle orbitals. But in answer to the question to whether it is 
possible to reconstruct a local operator if only its p-h matrix elements are given, in 1 10] it is 
remarked that it seems not possible except for a special case when the number of occupied 
particle states («) equals one. It is stressed by Holas et al. II24I that there is a positive answer 
for any n given in Sec. Ill of the work of Holas and Cinal 1251 . 

As a second focus arising from the work in LIOJ , Holas et al. 1241 refer to the differential 
equation satisfied by the density amplitude w(r) '' ^, involving the concept of the Pauli poten- 
tial. While conceptually satisfying, in f24l some reasons are set out why the implementation 
of such a radical Kohn-Sham scheme 1101 is presently hardly computationally competitive 
with that based on Slater-Kohn-Sham (SKS) orbitals. Being numerically demanding already 
for the static case, we expect the computational scheme discussed by Schirmer and Dreuw 
to be even more involved in the full time-dependent case. 



3 The differential virial theorem in time-dependent theory 

We turn next to an important result for DFT: namely the differential virial theorem (DVT). 
Following the earlier study of March and Young 1451 on the idempotent Dirac density matrix 
via its equation of motion, which led then in one-dimension to a result which they termed 
the differential form of the virial theorem, Holas and March ]46] some four decades later 
established the DVT in three dimensions with full account also taken of electron-electron 
interactions. Ref. ]46] has a fairly direct generalization to time-dependent theory 113l ll4lfT5] 



l40l and leads to the result (see also HI): 
d n{r,t) _ lv4„(^_y)^y^(^^j)^2V / dr'mir.r' ,t)Vu{r,r')+V\n(r,t)VV,^,{r,t)] 



dt^ 4 

(14) 

Here H2(r,r',f) is the pair density, M(r,r') is the electron-electron repulsion potential energy, 
while z(r, f ) is a vector field defined from the kinetic energy density tensor tafi (r, t) follow- 
ing Holas and March II46I . To be a little more specific, eqn(ll4b is, in essence a combination 
of the DVT and the continuity equation relating density n(r,f) and current density j(r,?), 
namely 

^ + Vj(r,0=0, (15) 

the latter being already invoked in the early work of March and Tosi [1] on TDDFT. What 
we want to stress is that eqnl ll4b can be employed at least in principle to construct the ap- 
plied (external) potential from the density n(r,f) plus the initial conditions, a matter already 
touched on in Sec.[T] Prerequisites for such a construction are accurate approximations for 
the kinetic energy density tensor and the correlated pair density. Efforts to obtain the latter 
in terms of first-order density matrices \A1\, should be very useful in this respect. 
Rewriting eqn(ll4li for the non-interacting SKS system, we have 

^-^ = -ivVr,f)+Vz,(r,f)+V[,7(r,0Vy(r,r)], (16) 

which requires only the knowledge of non-interacting vector field Zs(r,r) 1 14] to extract the 
Kohn-Sham potential V(r,?) from a given density n(r,r). This should be helpful in cases 
where direct inversion of the Kohn-Sham equations is impossible (cf. sub-section |2.1| i. 

In order to obtain expressions for the exchange-correlation potential as functional of 
the density, admittedly, the pair density n2(r,r',?) still enters in the subtraction of yext(r,r) 
from V (r, f ) . In Ref . 1 46 1 dealing with the time-independent problem, the correlated pair den- 
sity and kinetic energy density tensor were replaced by their non-interacting counterparts, 
which lead after a combination of the analogues of eqnl ll4b and eqn(|16t to an exchange- 
only approximation of the exchange-correlation potential beyond the Slater form. A similar 
proceeding is expected to succeed also in the present time-dependent case. As discussed 
by Qian and Sahni ll3lll4lfT5l in their derivation of eqnl ll4l l and ( I16l l. the combination of 
the above equations is also important from a more formal point of view, since it allows to 
disentangle several contributions to the exchange-correlation functional and interpret these 
in physical terms. 

4 Shortcomings of present calculations by TDDFT on charge-transfer excitations 

It is highly relevant to the foundations of TDDFT that current usage leads to substantial 
errors for charge-transfer excited states |48l|49ll50l . Usually, the excitation energies are 
severely underestimated. Furthermore, the potential energy curves of such charge transfer 
states do not display the known \/R dependence along a charge-separation coordinate R 1511 
I51531. 

A long-term solution to this problem may well lie in the use of time-dependent current 
density functional theory, which has recently been implemented |54| (see also |1 3,53])- 
It turns out, as discussed in |54 1 that a correct description of charge transfer excited states 
requires non-locality, and current density functional have this property. However, the high 



computational cost of such a current density approach raises doubts as to whether this route 
will be applicable to large molecules in the foreseeable future. 



5 Summary 

Serious criticisms of the foundations of TDDFT have recently been made flOl, the focus 
being on parts of the Runge-Gross work |8|. The example set out in sub-section 12. II and 
culminating in eqn(|7]l, shows in an admittedly simple time-dependent problem with an exact 
solution, that the time-dependent density can be correctly calculated from a one-body time- 
dependent potential in support of the mapping theorem |8|. 

We stress here again that the mapping theorem is not challenged in 1101 . The major 
point of the Schirmer-Dreuw study is rather the claim that the time-dependent Kohn-Sham 
approach has no predictive power due to the lack of a valid variational principle. How- 
ever, once the time-dependent density, discussed in some detail in sub-sections 12. 1 1 and 12.21 
above, is obtained by some other means, and the exact time-dependent exchange-correlation 
potential-functional would be known, the exact time-evolution of the electron density of the 
interacting system can be reproduced by the time-dependent Kohn-Sham equations. Since 
this manuscript was completed, a comment |56 J and reply |57 1 on the Schirmer-Dreuw study 
1101 have been published. 

Further matters discussed involve questions of the causality of the potential in TDDFT, 
first raised by Gross and co-workers I34II35I and subsequently discussed by Amusia and 
Shaginyan |36,37] and Harbola |40], and of the need to face the additional complications 
of current density theory for a specific class of excitations, namely charge-transfer excited 
states. This is because of the fundamental need for non-locality, which is correctly embedded 
in current density theory. Questions then arise as to the feasibility of application of such an 
approach to large molecules, because of the high cost. 
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